Magneto-transport Properties near the Superconductor-Insulator Transition in 2D 
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We analyze here the behavior near the 2D insulator-superconductor quantum critical point in the 
presence of a perpendicular magnetic field. We show that with increasing field H, the quantum disor- 
dered and quantum critical regimes, in which vortex degrees of freedom are suppressed, crossover to a 
new magnetically activated (MA) regime, where the correlation length ^ ~ 1/y/H. In this regime, we 
show that the conductivity decreases monotonically as opposed to the anticipated saturation predicted 
from hyperuniversality arguments. This discrepancy arises from the lack of commutativity of the fre- 
quency and temperature tending to zero limits of the conductivity. In the low-field regime such that 
<^ A, and in the absence of Ohmic dissipation, where A is a measure of the distance from the 
quantum critical point, the resistivity saturates to the Bose metal value found previously for Cooper 
pairs lacking phase coherence. 
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From the earliest experiments on superconduct- 

ing thin metal alloy films, a simple picture emerged 
for the magnetic-field tuned transition to the insulating 
state. Namely, the transition is continuous with a well- 
defined critical field. He, below which superconductivity 
occurs and above which an insulator obtains. At the 
critical field, the resistivity is independent of tempera- 
ture. Within the dirty boson model [Q in an applied 
magnetic field, these results follow naturally as a direct 
consequence of the inherent duality between Cooper pairs 
and vortices: above He, vortices condense while Cooper 
pairs remain localized but below He it is the Cooper 
pair condensation and vortex localization that affords 
the zero-resistance state. At criticality. Cooper pairs on 
the brink of forming a phase coherent state diffuse with 
the quantum of conductance ctq = Ae^/h for charge 2e 
bosons. However, more recent experiments have shown 
that this simple picture breaks down dramatically at and 
below He- In particular, a metallic phase HQ emerges 
below He with the zero-resistance state occurring only 
at a much smaller field Hm ^ He On phenomeno- 
logical grounds, it has been suggested that the metallic 
phase arises from the quantum tunneling of field-induced 
vortices in the presence of dissipation . 

As the metallic phase has been observed at zero-field 
pO| in the disorder-tuned insulator-superconductor tran- 
sition as well, we reinvestigated the conductivity of the 
standard quantum-disordered regime in zero field. We 
found that quite generally Cooper pairs lacking phase co- 
herence form a Bose metal phase This result arises 
from the fact that in the quantum-disordered regime 
the elementary excitations obey Boltzmann statistics. 
Hence, their concentration n is exponentially suppressed. 
As a consequence, so is the inverse scattering time 1/r 
that is determined by the probability that two quasipar- 
ticles collide. But because the conductivity is a product 
of the density and the scattering time, the exponentials 



cancel one another, thereby leading to a finite conductiv- 
ity as T ^ 0. In the current study, we focus on the fate 
of the Bose metal phase in a magnetic field. Near the 
quantum critical point where vortex dynamics are sub- 
dominant, the relevant physics can be accessed by the 
standard 0{M) Ginzburg-Landau (GL) theory In 
this work, we utilize this approach for a system near the 
1ST in the presence of a perpendicular magnetic field. 
We establish that at low temperature T and close to 
the quantum critical point, A ^ 1, the application of 
a field H drives the system to the high-field magnetically 
activated (MA) regime once Vi7 ^ A,T, thereby com- 
pounding the analysis of the transport properties in the 
quantum disordered regime. In this regime, the energy 
levels of the bosonic quasiparticle excitations are quan- 
tized and have a gap m ~ ^/H that is inversely propor- 
tional to the correlation length ^. Working within the 
large M approximation, we discuss the crossovers and 
limiting forms of the scaling functions for the correlation 
length ^ = ^{H, T, A) and the static conductivity a in the 
vicinity of the critical point. We present the qualitative 
form for the conductivity in the presence of dissipation 
based on the standard Kubo formula neglecting the ef- 
fects of frustration. The dc conductivity is regularized 
by inclusion of the weak coupling to an Ohmic heat bath 
rj as well by quasiparticle scattering 1 /r due to the quar- 
tic non-linearity in the GL action. Collectively, internal 
and external dissipation can be treated jointly with an 
effective coupling constant fj = rj + 1/r. In the absence 
of the bath and at a small constant magnetic field such 
that \/H <C a, the magnetic energy scale is irrelevant. 
Consequently, the Bose metal obtained previously in 
the quantum disordered regime is recovered. We show 
that in the MA regime, the conductivity cr is directly, 
rather than inversely, proportional to the total dissipa- 
tion 77. Further, in contrast to the low-field limit, the 
conductivity in the MA regime is found to be a decreas- 
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ing function of temperature and inversely related to the 
field. 

Our starting point is the M-component Ginzburg- 
Landau free energy |15|,p^, 
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in the presence of a magnetic field, where Ao(r) = 
{0, Hx, 0} is the vector potential due to the applied field, 
e* = 2e, and S is the T — unrenormalized distance to 
the quantum critical point. For a time- independent mag- 
netic field, it is convenient to expand the order parameter 
Tpa{ic,T) in the corresponding eigenfunctions that diago- 
nalize the Gaussian part of the free-energy density 
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where 5'„(a: — hpy/e*H) is the solution of the corre- 
sponding eigenequation, expressible in Hermite polyno- 
mials, and Ca(n,u;m,Py) are the expansion coefhcients. 
The related eigenvalues are given by LOuin + 1/2) with 
ljh = 2e*Ha/h — Airf /a, where a is the lattice constant 
and / = <i>/$o is the magnetic frustration. Henceforth, 
we rescale H such that uh H . In the large M (mean- 
field) treatment the bare T = distance to the critical 
point 5 is renormalized, and the equation for the quasi- 
particle excitation gap m, 
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must be solved to obtain the effective Gaussian action. 
Performing the summation over the Matsubara frequen- 
cies exactly and employing the Poisson formula for the 
sum over 'Landau' levels n, we obtain (assuming m <C 1) 
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FIG. 1. Phase diagram for the field-tuned transition in the 
presence of quantum fluctuations, S. The zero-temperature 
quantum critical point studied here corresponds to 5 — 5c. 
Quantum disordered (QD), quantum critical (QC) and mag- 
netically activated (MA) regimes correspond to the limiting 
cases A > \fH,T, T > \/h,A and \/!ff > T, A, respec- 
tively. 



Three limiting regimes follow from Eq. (^. 1) If 
T, ViJ <C A, then m — A and we are in the quantum 
disordered regime (QD). This regime has been studied 
extensively recently, the key surprise being the onset of a 
Bose metal phase once the Cooper pairs loose their 
phase coherence. 2) The regime \/H, A <C T is quantum 
critical (QC) in which m = OT (6 = 2ln{{V5 + l)/2)). 
Extensive studies of this regime exist in the literature 
]l7t . In both of these regimes, the magnetic field effects 
are subdominant, leaving the conclusions of the previous 
analyses [pT|JT7| in tact. 3) Finally, when y/H > T, A, 
the gap m = AVh emerges, where A — 0.418945 is the 
solution to the equation 
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and A = (At:/U)6 + A (A is the upper momentum cut- 
off) is the renormalized distance to the zero-temperature 
critical point in the absence of the field. 



In this regime, which we refer to as magnetically acti- 
vated (MA), the correlation length, m~^, exhibits univer- 
sal behavior and is determined (in the leading approxi- 
mation) solely by the magnetic field H. The energies 
of bosonic quasiparticles, which represent the low-energy 
excitations close to the transition point, are quantized 
and occupy the degenerate 'Landau levels' with energies 
e„ = V Hn + m?. The 1/M correction is expected to 
modify the universal value of A. Its calculation is pro- 
hibitively difficult because the coefficient in the quartic 
term of the GL action is a complicated function of the 
energy levels e„. All of these regimes imply that the ele- 
mentary excitations are spin-wave like. This takes place 
only if A > 0. In the opposite case, close enough to the 
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finite T, finite H transition, transport is determined by 
the motion of vortices as illustrated schematically in Fig. 
(|l|). The finite-size scaling considerations suggest that 
the correlation length near the transition is a universal 
function of two arguments: £, = T'^ f{A/T, VH /T) and 
assumes the appropriate limiting forms if one of the pa- 
rameters tends to zero. The M = oo form of m = 
as a function of \JillT obtained from the numerical so- 
lution of Eq. for different values of A/T is presented 
on Fig. (H). As is evident, a crossover occurs to Vi? 
dependence as the field increases. 

The linear conductivity (per flavor) near the 1ST can 
be obtained from the GL action with the help of the Kubo 
formula 
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While calculating the averages in the partition function 
we expand ip using Eq. which yields for the longitu- 
dinal conductivity |18|]: 
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The experimentally observable finite temperature con- 
ductivity must be calculated once dissipative mecha- 
nisms are taken into account | p5[ . Weak dissipation, 
fj = rj + l/r, is accounted for by introducing the term 
?7|cli„i||V-'(i", into the free-energy density Its 
inclusion changes the analytical properties of the Green 
functions in Eq. (||) . Performing the analytical con- 
tinuation, one finds that the static conductivity is given 

by 
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The analysis of Eq. (Q) simplifies if the dissipation is 
weak compared to the quasiparticle energies, implying 
rj <^ m. Asymptotic analytical expressions are readily 
obtained for m ^ T, which define QD and MA regimes. 
In these cases, there are two competitive contributions 
to the integral over x. One comes from the minima of 
the denominator at x^ « {e"^ — ff /2)/AT'^ and x'^j^^ k. 
(4+1 - ?7V2)/4r2 and the second fr om the region near 
X = Q. Performing the approximate integration over x, 
and analyzing subsequently the sums over n, we arrive 
at the following limiting cases. 



In the QD regime, m = A ^ \/H, and the conductiv- 
ity is given by cr = cr'-^-' -I- cr^'^\ where 
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AT < i? < A2 



If we are in the MA regime with \/ H ^ m, 
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In the formula above, 



n + l 



{n + 1 + A^)^n + A^y 



23.913332 



is a rather large numerical prefactor, arising from the 
M = 00 limit. These expressions describe qualitatively 
the crossover from the QD to MA regime, provided the 
quasiparticles are well defined- that is, dissipation is suf- 
ficiently weak. As the field increases, the conductivity 
becomes directly proportional to dissipation, as dictated 
by Eq. (|^) . This can be attributed to the splitting of the 
minima in the denominator of the general formula, Eq. 

There are two sources of dissipation for which we must 
account. Internal dissipation arises from mutual scatter- 
ing of quasiparticles, leading to a finite quasiparticle scat- 
tering rate l/r. External Ohmic dissipation 77 may come 
from the resistive shunting of superconducting grains in 
a system of fabricated Josephson-junction arrays (JJA). 
The coefficient r/ is then inversely proportional to the re- 
sistivity of a shunting resistor, Rg [po| , which may be 
itself a function of temperature and magnetic field. In 
addition to Ohmic dissipation, disorder can also serve as 
a dissipative channel as has been shown for static disorder 
in the presence of d-wave pairing ||2^ . In the QC regime, 
interactions alone lead to a large inverse scattering time, 
1/t ^ T ^ m. Hence, quasiparticles are poorly defined 
jisf , and small e = 3 — d or 1/M expansions are neces- 
sary to determine the conductivity in a controlled way 
]l2t . The implementation of this procedure leads to the 
universal value for the conductivity close to aq — Ae'^/h 
]l5[ . The weak external dissipation <ti m T is clearly 
subdominant in this regime, except at low temperatures 
]l6[ . On the contrary, in the QD regime, l/r ^ Te"^^'^, 
(up to logarithmic corrections), as a result of the Boltz- 
mann statistics of the quasiparticles. In the absence of 
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Ohmic dissipation, as follows from Eq. (||)- metal- 
lic, rather than insulating behavior of the conductivity 
obtains provided that -s/ff <C A |11 . The persistence of 



the Bose metal phase for a finite range of magnetic field 
is possibly relevant to the recent experiments of Mason 
and Kapitulnik 0. 





FIG. 2. The dependence of m/T and the corresponding 
conductivity a (in units of aq = 4e^//i) on \fli jT for 
three values of the distance to the quantum critical point: 
A/T = (solid line), A/T = 1.0 (dashed line) and A/T = 3.0 
(dashed-dotted line). Increasing A results in a crossover from 
the quantum critical (A <^ T) to the quantum disordered 
(A > T) regime. 



In the MA regime, quasiparticles also obey Boltzmann 
statistics. Based on dimensionality arguments similar to 
those valid for the QD regime, we suggest that in the 
MA regime, 1/r ^ Te^-^^^'^ , where the constant of 
proportionality is determined by a numerical coefficient 
that is of the first order in 1/M. The calculation of the 
corresponding lifetime is complicated because of the dis- 
crete character of the quasiparticle spectrum. However, 
as follows from Eq. (^), the conductivity becomes di- 



rectly, rather than inversely, proportional to dissipation 
for sufficiently large H. This means that inclusion of 
external dissipation i] and determines only how fast the 
conductivity decreases with increasing field. To illustrate 
this, we resort to a semi-phenomenological calculation in 
which we set 1/r = 2.3300502 • Te"™/'^ and 77 = O.IT. 
The interpolation formula for 1 /r is chosen to satisfy the 
behavior in all regimes near the critical point. The nu- 
merical prefactor in 1 /r is taken in such a way, that for 
?7 = 0, cr ~ (Tq in the QC regime. The conductivity was 
then calculated numerically as a function of \/H/T for 
different values of A/T, using Eqs. (^ and The re- 
sults that are depicted in Fig. (0) support rather rapid 
decrease of a with increasing \/H/T but their depen- 
dence on rj/T is relatively weak. We conclude, then, 
based on the mean-field treatment, that the conductivity 
decreases when the crossover to the MA regime occurs 
even if 77 = 0. This behavior is in agreement with ex- 
periments on the magnetic field-tuned 1ST in disordered 
indium-oxide thin films (l). However, these results are in 
disagreement with the suggestion of Kim and Weichman 
that the conductivity saturates with increasing field . 
This discrepancy can be explained by the fact that for 
?7 = in the general universal scaling function, 

a = (4eV/i)S(cj/r, A7T, Vi//A), (11) 

applicable for the conductivity near the d = 2 1ST, the 
limits T = Q, u! ^ and a; = 0, T — > do not com- 
mute The former is the coUisionless limit, appar- 
ently considered in Ref. (14), in which the conductivity 
is likely to experience a crossover between two universal 
values depending on the ratio ^/H / A |l8j. However, it is 
the collision-dominated limit that is experimentally rele- 
vant and hence the main subject of our treatment. It is 
straightforward to verify that Eq. (Q) satisfies the general 
scaling relation ([ll| ) with uj — and z — 1, provided that 
77 = and the inverse scattering time 1/r obeys itself 
the scaling form 1/r = Tf{A'-/T, VH/A) in the critical 
region. 

The results above were obtained for low enough mag- 
netic fields, ^/H ~ / ^ 1 and hence, magnetic frustra- 
tion can be neglected. In translationally periodic systems 
such as J JA, magnetic frustration plays an important role 
giving rise to peaks in the conductivity for integer values 
of frustration / |2^. Moreover, the peaks correspond- 
ing to the rational values / = 1 1 1 1 f • • • may also be 
observed |Q. However, disorder disrupts the lattice pe- 
riodicity smearing the peaks in the conductivity giving 
rise to a disappearance of frustration effects. This situa- 
tion takes place in the disorder-tuned 1ST in thin films. 
Hence, we expect, that as long as Cooper pairs remain 
intact and the transition is described by Eq.(Q), our cal- 
culations provide a qualitative picture for the behavior 
of the magneto-conductivity near the point of nominal 
the 1ST in thin films including the observation of a Bose 
metal phase. 



4 



This work was funded by the DMR of the NSF grant 
No. DMR98-96134. 



A. F. Hebard and M. A. Paalanen, Phys. Rev. Lett. 65, 
927 (1990). 

A. Yazdani and A. Kapitulnik, Phys. Rev. Lett. 74, 3037 
(1995). 

For a review see A. M. Goldman and N. Markovic (1998) 

Physics Today 51 (11) 39-44. 

M. P. A. Fisher, Phys. Rev. Lett. 65, 923 (1990). 

Y. Liu and A. M. Goldman, Mod. Phys. Lett. 8, 277 

(1994). 

N. Mason and A. Kapitulnik, Phys. Rev. Lett. 82, 5341 
(1999). 

D. Ephron, A. Yazdani, A. Kapitulnik, and M. R. 
Beasley, Phys. Rev. Lett. 76, 1529 (1996). 



N. Mason and A. Kapitulnik, |gond-mat/0006138 . 
E. Shimshoni, A. Auerbach, and A. Kapitulnik, Phys. 
Rev. Lett. 80, 3352 (1998). 

H. M. Yaeger, et. al., Phys. R ev. B 40, 182, (198 9). 
D. Dalidovich and P. Phillips, |cond-mat/0005l"l£ , to ap- 
pear in PRB, 2001. 

Subir Sachdev, Quantum phase transitions (Cambrige 

University Press, New York, 1999). 

Min-Chul Cha, et. al. Phys. Rev. B 44, 6883, (1991). 

K. Kim and P. B. Weichman, Phys. Rev. B 43, 13583 

(1991). 

K. Damie and S. Sachdev, Phys. Rev. B 56, 8714 (1997). 
D. Dalidovich and P. Phillips, Phys. Rev. Lett. 84, 737, 
(2000). 

A. Chubukov, S. Sachdev and J. Ye, Phys. Rev. B 49, 
11919 (1994). 

A. van Otterlo, K.-H Wagenblast, R. Fazio, and G. 
Schon, Phys. Rev. B 48, 3316 (1993). 
A. O. Caldeira and A. J. Leggett, Phys. Rev. Lett. P46, 
211 (1981). 

K. Wagenblast, A. van Otterlo, G. Schon, and G. Zi- 
manyi, Phys. Rev. Lett. 79, 2730 (1997). 

H. M. Jaeger, D. B. Haviland, B. G. Orr, and A. M. 
Goldman, Phys. Rev. B 40, 182 (1989). 

I. F. Herbut, Phys. Rev. Lett. 85, 1532, (2000). 

R. A. Webb, et. al. Phys. Rev. Lett. 51, 690, (1983) 
H.S.J, van der Zant, et. al. Phys. Rev. B 54, 10081, 
(1996). 



5 



